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We investigate some aspects of the renormalization group flow of gravity in the presence 
of fermions, which have remained somewhat puzzling so far. The first is the sign of the 
fermionic contribution to the running of Newton's constant, which depends on details of the 
cutoff. We argue that only one of the previously used schemes correctly implements the 
cutoff on eigenvalues of the Dirac operator, and it acts in the sense of screening Newton's 
constant. We also show that Kahler fermions give the same contribution to the running 
of the cosmological and Newton constant as four Dirac spinors. We then calculate the 
graviton contributions to the beta functions by imposing the cutoffs on the irreducible spin 
components of the tetrad. In this way we can probe the gauge dependence of the off-shell 
flow. The results resemble closely those of the metric formalism, except for an increased 
scheme- and (off shell) gauge-dependence. 



I. CUTOFF SCHEMES 



The functional renormalization group equation (FRGE) is a convenient method to calculate 

n □, 

the quantum effective action that is not restricted to renormalizable field theories [1|, • Instead of 
computing directly the effective action, one computes its derivative with respect to an external mass 
scale k called "the cutoff". This "beta functional" is automatically free of UV and IR divergences. 
The advantage of this procedure is that it provides a definition of the beta functions of a theory 
that does not refer to any UV regularization. Then, taking some (generally local) action as starting 
point at an UV scale A, one can obtain the effective action by integrating the FRGE from A to zero. 
The divergences of the theory reappear when one tries to shift A to infinity and can be analysed 
by integrating the FRGE in the direction of increasing k. 

For our purposes the most useful implementation of this idea is the Wetterich equation 
It gives the ^-dependence of a functional called the Effective Average Action (EAA), which is a 
modified effective action where the propagation of the low momentum modes is suppressed. Here 
"low momentum" generally means "momenta below a given scale k". The suppression is described, 
in flat space, by some function Rh(q 2 ), where q 2 is the squared momentum. To be admissible, the 
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cutoff function must satisfy the following basic requirements: (i) it must be continuous and 
monotonically decreasing in q 2 and k, (ii) it must go rapidly to zero for q 2 > k 2 , (iii) it must 
tend to a positive (possibly infinite) value for q 2 — > and (iv) it must tend to zero for k — > 0. 
It is unavoidable that details of the flow depend on the shape of the cutoff function, even though 
the effective action, which is obtained by letting k — > 0, is not. This is referred to as "scheme 
dependence" and is akin to the renormalization scheme dependence in perturbation theory. 

When one considers anything more than a single scalar field, in addition to the shape of the 
function Rk(z), one encounters further choices. We will be concerned mainly with gravity, which 
offers a particularly complex scenario. The application of the FRGE to gravity has been pioneered in 
Q|, followed by Q|. Already in these early works the implementation of the cutoff followed different 
methods. In [3j the cutoff was imposed on the tracefree and trace parts of the metric fluctuation (a 
decomposition that is purely algebraic), whereas in (4] the cutoff was imposed separately on each 
irreducible spin component (a finer decomposition that involves differential conditions). Following 
[5(1, we will call the former a cutoff of type "a" and the latter of type "b". When gravity is involved, 
the kinetic operators governing the propagation of quantum fields typically have the structure 
—S/ 2 +aR+bE, where V is the Riemannian covariant derivative, R stands here for some combination 
of curvature terms (Riemann, Ricci and scalar curvatures, with indices arranged in various ways) 
and E a term that may involve other background fields and possibly also couplings. Since Fourier 
analysis is not available in general curved spaces, in defining the cutoff one has to specify an operator 
whose eigenvalues form a basis in field space and in general play the role that plane waves play 
in flat space. One takes Rk to be a function of this operator (or rather, its eigenvalues). In [5fl 
the following terminology has been introduced: a cutoff is of type I if R^ is a function of —V 2 , of 
type II if it is a function of —V 2 + aR and of type III if it is a function of —V 2 + aR + bE. It is 
clear that there are infinitely many more possible cutoff variants and that this classification is very 
incomplete, but it covers the most commonly used cutoffs. 

The RG flow in the so-called Einstein-Hilbert truncation, where one retains in the action terms 
at most linear in curvature, has been studied in (g - 10 j , with en extended ghost sector in 



and u sing the Vilkoviski-de Witt formalism in 
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13] 



14l | . Higher derivative terms have been analysed in 



20fl. The contribution of matter fields has been taken into account in 



I cutoff and in 
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24l | using a Type 



25| usin g a 



has been studied in 



26r|29fl. 



type II cutoff. Conversely, the effect of gravity on fermionic interactions 



Physical observables, which are related to the full effective action at k = 0, will be independent 
of the cutoff choice. Furthermore, some terms in the beta functional, namely those that refer 
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to dimensionless couplings, can actually be shown to be scheme-independent. It is also well- 
understood that if the system has a nontrivial fixed point, its position is not universal. In general, 
the scheme-dependence of the flow should therefore not be a cause of major concern 30]. 

In the literature on the renormalization group for gravity, there is however one result where the 
scheme dependence seems to be particularly nasty: it is the contribution of fermion fields to the 
beta function of Newton's constant. This contribution has been computed for example in section 
III of 



ij. If there are n d Dirae lields il is 
if one uses a type I cutoff and 

- / <*- ^ G - 1) ft (^) 4 * = (^) / ** * * <*) 

if one uses a type II cutoff. Here Q n {W) = jJU f °° dzz n ~ 1 W{z). If we use the optimized cutoff fl 

R k (z) = (k 2 - z)6{k 2 - z) , (3) 

we get Qi ) = 2k 2 , Q2 (^ 9i p2 k ^ = k 2 and therefore the contribution of fermions to the relevant 

term of the beta functional is 

n D 



96vr' 



j d^x^/gR for the type I cutoff (4) 

2 / d 4 XyfgR for the type II cutoff. (5) 

n 

The latter agrees with an earlier independent calculation of the renormalization of G in [3l| : the 
former differs not just in the value of the coefficient but even in the sign. 1 To check that this is 
not just a quirk of the optimized cutoff, consider an exponential cutoff 

ze ~a(z/k 2 ) b 

R ^ = I _ e -a( Z /k*)» ■ ( 6 ) 

In order to guarantee that condition (iii) is satisfied, one has to assume b > 1. For example if b = 1, 
Qi (j^pj^j = Tr 2 k 2 /3a and Q2 ^^ fc J = 2k 2 /a and therefore we encounter the same problem, 
independently of the parameter a. One can see that the same is true also for b > 1. 

This sign ambiguity is puzzling. Since the beta function of G vanishes for G = 0, the sign of 
G can never change in the course of the RG flow (we disregard here as pathological the case when 



1 This was noted a while ago by Calmet et al. while they were working on (33 |. R.P. wishes to thank D. Reeb for 
correspondence on this point. 
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the inverse of G passes through zero). Therefore, if the physical RG trajectory approaches a fixed 
point in the UV, the sign of Newton's constant at low energy will be the same as that of Newton's 
constant at the fixed point. In a model where gravity is induced by minimally coupled fermions, 
the sign of Newton's constant would depend on the scheme. 

One has to be careful in drawing physical conclusions from these calculations: the relation 
between the coupling G appearing in and the physical Newton's constant that is measured in 
the lab may not be as simple as it seems. The functional that obeys an exact RG equation is 
not a simple functional of a single metric but rather a functional of a background metric and (the 
expectation value of) the fluctuation of the metric. This functional is not invariant under shifts 
of the background and fluctuation that keep the sum constant, and therefore there may be several 
couplings that one could legitimately call "Newton's constant" (for a discussion see 33j|). The 



one that we are discussing here is the one that multiplies the Hilbert action constructed from the 
background metric only, and it is not obvious that the other ones would behave in the same way. In 
order to make a completely well-defined statement one should really calculate a physical observable 
and in such a calculation all ambiguities should disappear. It is therefore possible, in principle, that 
even the sign difference between the two calculations discussed above may be resolved when one 
considers physical observables. In this paper we will argue for a different and simpler solution of 
the issue, namely that only the type II cutoff gives a result with the physically correct sign. 

As an aside, we will also calculate the contribution to the beta functions of the cosmological con- 
stant and Newton constant due to Euclidean Kahler fermions in four dimensions. Kahler fermions 
are a way of describing fermion fields in terms of Grassmann algebra elements, instead of spinors 
and has the merit that it does not require the use of tetrad fields. We find that one Kahler fermion 
gives exactly the same contribution as four Dirac spinors, and that the same sign issue is present. 

When gravity is coupled to fermions an additional question arises regarding the field carrying 
the gravitational degrees of freedom. In jj, '. 
added to that of matter fields, as well as in 



3: 



22 L where the contribution of graviton loops was 



26H28| where scalars and fermions were interacting 
with gravity, the carrier of the gravitational degrees of freedom was the metric. It has been pointed 
out in 



34| that in the presence of fermions it may be more natural to use the tetrad. Even if one 
chooses a Lorentz gauge where the antisymmetric part of the tetrad fluctuation is suppressed, the 
two calculations are not the same, because one has to work off shell and the Hessian in the tetrad 
formalism contains some additional terms proportional to the equations of motion. The calculations 

n 

in [34], which used a type la cutoff in the a = 1 gauge, show that a fixed point is still present, but 
is less stable than in the metric formulation. We will compute the gravitational contributions in 
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the tetrad formalism, using type b cutoffs (both lb and lib); this allows us to analyse the off-shell 
a-dependence. The results will be found to be somewhat closer to those of the metric formalism 
than those found in [34] (which we have independently verified), but they still present a stronger 
dependence on the scheme and on the gauge parameter. 

II. KAHLER FERMIONS 

Before discussing in detail the issue of cutoff scheme with ordinary spinor fermions, we would 
like to point out that precisely the same problem arises also in another representation of fermionic 
matter. This section is not strictly necessary for the following discussion so readers that are mainly 
interested in the solution of the puzzle presented above may go directly to section III. 

As is well-known, in any dimension the Grassmann and the Clifford algebras are isomorphic as 
vector spaces. This is the basis for a representation of fermion fields as inhomogeneous differential 



forms [35l437l|. Such fields are called Kahler fermions. In this representation the analogue of the 
Dirac operator is the first order operator d + 5, where d is the exterior derivative and 5 is its 
adjoint. Note that the definition does not require the use of a tetrad. We would like to compare the 
contribution of Kahler fermions to the gravitational beta functions to the one of ordinary spinor 
fermions. In particular we would like to see whether the same sign issue arises when different cutoff 
types are used. Since the details of the calculation are strongly dimension-dependent, we shall 
restrict our attention to the case d = 4. 

An inhomogeneous complex differential form $ can be expanded as 

<I> = if(x)+ip^(x)dx^+—^^ u (x)dx^f\dx u +—^^ up (x)dx^f\dx u f\dx p +—^ 
.2 . 3. 3. 

We can map the 3- and 4-form via Hodge duality into a 1- and 0-form respectively. The field <3? 
thus describes a scalar, a pseudoscalar, a vector, a pseudovector and an antisymmetric tensor, for a 
total of 16 complex components. This is an early sign of the fact that one Kahler field is equivalent 
to four Dirac fields. 

The square of the Kahler operator is the Laplacian on forms: A = (d + 5) 2 = dd + 5d. On forms 
of degree 0, 1 and 2 it is given explicitly by 

A<°> = -V 2 (7) 
A (i)/i = _v 2 ^ + R\ (8) 

A (2) 7CT = _ V 2 1q/?7ct + R j Sf}l , _ R ^a _ 2R J * (g) 
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In order to read off the beta functions of the cosmological constant and Newton's constant it enough 
to consider a spherical (Euclidean de Sitter) background, with curvature tensor 
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w " d(d 
Then the operators defined above reduce to 

AWf = l-v 2 + -p ) sn 

. (2)7(7 



P 



[IV 



d 



9fiiyR 



A 



a/3 



1 



1 



-V* + -P I 1, 



70" 



(10) 

(11) 

(12) 



The ERGE for a Kahler fermion with a type II cutoff is 
dT k 



dt 



2-Tr 
2 



(0) ^ P fc (A(0)) J " 2 ^ ^(Ad)) 



-Tr 



(2) 



P fc (AP)) 



1 1 



2 (4 



For a type I cutoff we find instead: 



-4Q 2 



~pT 



dt 



1 f d t R k (-V 2 ) 



2-^Tr (1) 



5 t i? fe (-V 2 ) 



(13) 



:4-i- 1 



2 (4vr) 5 



d^Xyfg 



-4Q 2 



dtRk 



P fc (-V 2 ) + f 



Pi 



-Tr 



Q2 



/ c\P fc (-V 2 ) 

(2) ^ fe (-V2) + f 

d t R k 

p2 



Evaluating the Q-functionals for the cutoff ([3]) we get: 



Type II: 
Type I: 



dt 
dT k 



4- 



1 1 



2 (4tt) 
4.1- 1 



j d A x^/g 



-Ak 4 + -Rk 2 



(14) 

(15) 
(16) 



dt ' 2(4vr) 2 

In both cases the effect of one Kahler fermion is seen to match exactly the result of four spinors 
(no = 4). This should not induce one to believe that spinors and Kahler fermions are equivalent: 
in fact, their contributions to the curvature squared terms are quite different. Nevertheless, the 
puzzling sign issue of the P-term that afflicts spinor fermions is present in this case too. 



III. CUTOFF CHOICE FOR FERMIONS 



We now return to ordinary Dirac spinor fields and we reexamine in more detail their contribution 
to the gravitational effective action and beta functions. For the sake of generality we will now work 
in arbitrary dimension d. The standard way of defining the effective action for a fermion field that 
is minimally coupled to gravity is to exploit the properties of the logarithm and write 



r 



-Trlog(|0|) 



iTrlog(0 2 ) 



1, 



Trlog -V 2 + 



P 



(17) 
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The corresponding EAA can then be defined as 

Tfc = -^Trlog 



o R 
-V 2 + - + R k 



(18) 



In the definition of this functional one encounters the same ambiguities that we have mentioned 
earlier for bosonic systems. In addition to the shape of the cutoff function R^, one seems to also 
have the freedom of choosing the argument of this function to be either —V 2 (type I cutoff) or 



V 2 + -? (type II cutoff). The former choice has been made in 



[4, 



21 



2J, the latter in |25|. Taking 



the £-derivative (where t = log(k/ko)) and defining Pk(z) = z + R]~(z), one has 

for a type I cutoff 



dt 



dt 



1 W-v 2 ) 

2 r ^(-V 2 ) + | 

l T ^i4(-V 2 + f) 
2 " P fc (-V 2 + f) 



for a type II cutoff. 



(19) 
(20) 



The first few terms in the curvature expansion of these traces can be evaluated, for any background, 
using heat kernel methods. However, for a spherical background, the spectrum of the Dirac operator 
is known explicitly and the same traces can also be computed directly as spectral sums. Comparison 
of these calculations indicates that only the type II cutoff correctly reflects the cutoff on eigenvalues 
of the Dirac operator. 

A. Heat kernel evaluation 



With a type I cutoff, the trace (|19|) giving contribution of njj Dirac spinors to the FRGE is 



dt 



n D 



2 [d/2] 



2 (47r) d / 2 



d t Rk 
Pk 



p, 



k / 



• (21) 



Here 2^- d / 2 ^ (where [x] is the integer part of x) is the dimension of the representation. The first term 
proportional to R is proportional to the heat kernel coefficient &2(~ V 2 ), and the second comes from 
the expansion of the denominator in (j 1 9 H . Evaluating the Q-functionals with the cutoff ([3]) we 
obtain 

2 {d/2] 



dl\ l_ 

dt r (| + 1) (47r) d / 2 ' 



n-D d Xy/g 



k a + 



d-3. 
~12 



k d ~ 2 R 



(22) 



Evaluating (|20p with the same techniques yields 



d£k 
dt 



n D 



2 [d/2] 

T (47r) d / 2 



+ ... 



(23) 



s 



where the term proportional to R comes entirely from the heat kernel coefficient 62 (—V 2 + -j). 
Evaluating with the cutoff ([3]) we obtain 



dt 



24 



(24) 



In d = 4 this yields the results quoted in section I. We see that the sign issue is present in any 
dimension d > 3. 



B. Spectral sums on S d 



The heat kernel calculation of the preceding subsection can be done in an arbitrary background. 
On the other hand, in the case of the spherical background we know explicitly the spectrum of the 
Dirac operator: the eigenvalues and multiplicities are 



A 



R 



d(d-l) \2 



d 



+ n 



m. 




n = Q,l,... 



(25) 



With this information one can compute the trace of any function of the Dirac operator as Tr /(0) = 
Yln^o m nf(^n)- We will now evaluate the r.h.s. of the FRGE by imposing a cutoff on the eigenvalues 
of the Dirac operator. The EAA can be defined directly in terms of the Dirac operator as 



r fc = -trio g (|#| + /#(|#|)) , 



(26) 



where the cutoff R® has to be a function of the modulus of the Dirac operator, since we want to 
suppress the modes depending on the wavelength of the corresponding eigenfunctions. This is also 
needed for reasons of convergence. Since the operator is first order, the conditions that R® has 
to satisfy are similar to (i)-(iv) of section I, except for the replacement of q 2 and k 2 by A n and k 
respectively. For the explicit evaluation, we will use the optimized profile 



R°(z) = (k — z)6{k - z) , [z > 0) . 



Then we have 



Tr 



p?{\lp\) 



E 



dtR£(\\n\) 



in, 



(27) 



(28) 



>niAni) ± „ 

The sum can be computed using the Euler-Maclaurin formula. Details are given in Appendix I. 
The result is 

1 



dt 



-Tr 



djRl 



r (I + 1) (4tt) 



(•(,/) ( k d -^k d ~ 2 R + 0(R 2 ] 



(29) 
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where V(d) is the volume of the tf-sphere. This agrees exactly with (|24p . which was obtained with 
the type II cutoff. (We have checked that the agreement extends also to the next order in the 
curvature expansion.) 

C. Discussion 

Note that computing the r.h.s. of the FRGE with a spectral sum is a much more direct procedure, 
since it avoids going through the square root of the square of the Dirac operator, and also avoids 
having to use the Laplace transform and the heat kernel. It is therefore also a more reliable 
procedure when there are ambiguities. The agreement of the spectral sum with the type II-heat 
kernel calculation is a useful consistency check and suggests that the latter gives the correct result 
whereas the type I cutoff does not. 

If so, there remain to understand why the type I cutoff should not be admissible in this case. 
One can get some hint by thinking of what this cutoff does in terms of eigenvalues of the Dirac 
operator. We begin by noting that (|26p can be rewritten as follows: 2 

T fc = -Itrlog (\lp\ + R£ (\lp\)) 2 = -itrlog (-V 2 + | + 2|0|i#(|#|) + Rg (|0|) 2 ) . (30) 

One can compare this with (118p . Note that the cutoff Rk in that formula could be a function of 
different operators which, on a sphere, differ by a constant shift. For the present purposes it is 
convenient to think of it as a function of ifi 1 '. Calling R^ this function and calling z = \Ip\, we have 

R k (z 2 ) = 2zR°(z) + R°(z) 2 . (31) 

We can solve this relation to get 

Rj?(z) = -z + Jz 2 + R k (z 2 ) , (32) 

so for any cutoff imposed at the level of (I18p one can reverse-engineer an effective cutoff to be 
imposed at the level of (|26p that will give the same result. 

In general, this cutoff may fail to satisfy the required conditions, in particular condition (iv). For 
a type II cutoff, Rk in (I18P is a function of z 2 , so Rk{z 2 ) = Rk{z 2 ). This implies that Ro(z 2 ) = 
and thus also Rq(z) = for all z > 0. For a type I cutoff, on the other hand, this may not be the 
case, as we will show in the following examples. 

2 This is a formal relation because the functional Tk is ill-defined, but one can write a corresponding relation for 
9tFfc, with the same result. 
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Consider first the optimized cutoffs. In the type II case one has Rk(z 2 ) = (k 2 — z 2 )9(k 2 — z 2 ) 
and one finds that in this case the corresponding cutoff R^{ z ) given by (|3"!Zj) is also optimized, and 
precisely of the form (|27p . This is a way of understanding why the two calculations give the same 
result. In the case of a type I cutoff, we have instead Rk{z 2 ) = (k 2 — z 2 + R/4)9(k 2 — z 2 + jR/4), 
whence one derives 



m*)=\ \k 2 + i~z 0[\k 2 + i-z . (33) 



This does not tend uniformly to zero when k — > 0. In the case of an exponential type II cutoff with 
Rk = Rk given by ([6]), we have 

R k D (z) = -z + Z , (34) 

yj\ _ e -a(z 2 /k 2 ) b 

which has all the desired properties. On the other hand for an exponential type I cutoff with R^ 
given by 




Rg(z) = -z + 



( R\ e-< z2 ~T )> 26 

f + V 2 -V ^7«^ - (35) 



For b odd, and in particular for the most natural case b = 1, this function does not tend uniformly 
to zero when k — > and therefore condition (iv) is not satisfied. 3 

These arguments lend support to the view that only the type II cutoff gives the physically correct 
result. Of course not all results obtained from the type I cutoff have to be wrong, for example the 
leading term (renormalizing the cosmological constant) and, in d = 4, the curvature squared terms, 
are the same using the two cutoffs. These however are just the "universal" quantities that do not 
depend on the choice of the cutoff. We believe that for the generic dimensionful, non-universal 
quantities, the results obtained via a type I cutoff should not be trusted. 



IV. TETRAD GRAVITY 



We will now compute the graviton contribution to the running of Newton's constant and cos- 
mological constant, in d-dimensions, when the tetrad is used as a field variable. This has been 

n 

discussed recently in [34] using a type la cutoff, i.e. a cutoff that depends on —V that is added 
to the full gravitational Hessian. In order to have a manageable, minimal Laplacian-type operator, 

3 Note that in the limit k — > the function is non zero only for z < y/ R/A. Since the smallest eigenvalue of the 
Dirac operator is ^/R/3, it remains true that limfc_>o I\- = T. 
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this requires that the gravitational gauge-fixing parameter be fixed to a = 1. We will use instead 
cutoffs of type b, meaning that the graviton is first decomposed into irreducible components of spin 
2, 1 and 0, and the cutoff is imposed separately one each component. This permits the discussion 
of general diffeomorphism gauges. We will use both type lb and type lib cutoffs. 



A. Hessian and gauge fixing 

The ansatz we make for the effective average action is the standard Einstein-Hilbert truncation 

T k [e,e] = T EH [e,e] +lf F [e,e] = J d d x dete (#(<?(e)) - 2A fc ) + lf F [ e ,e] (36) 

where we have indicated the ^-dependence of the couplings and T^ F is the gauge-fixing term, 
to be specified below. In tetrad formulation the metric is represented in terms of vielbeins as 
Qfiu = e a n,e b vVab- If we decompose g^ = g^ + h^ u and e a ^ = e a ^ + e a ^ we have the relation 

h,u, = 2e^ v) + £(/e„)p (37) 

where Latin indices on e have been transformed to Greek ones by contraction with e. Now sub- 
stituting this formula in the Taylor expansion of the metric in terms of metric fluctuations we 
find: 

/5T EH f 5F EH ICC 5T EH 

%7 2£ - + c / TST*"" + 2 / / « ^ + (38) 

In the third term we have introduced by hand a factor < £ < 1 that interpolates continuously 
between the pure metric formalism (£ = 0, h^ v = 2e^ v \) and the tetrad formalism (£ = 1). We 
see that the part of the action quadratic in e differs from the one in the metric formalism by terms 
proportional to the equations of motion. Since in the derivation of the beta functions it is essential 
to work off shell, we cannot ignore these terms. 

The gauge fixing terms for diffeomorphisms and local Lorentz transformations are 

rf F K e] = ^ j & A x^r v F»F v + ^- J d d x^ G ab G ab (39) 

For diffeomorphisms we choose the condition 

F " = 7K3(^-f^ S ) • (40) 

while for the internal O(d) transformation we choose a symmetric vielbein G ab = e^- ab h We will 
choose oil = in order to simplify the computation. This correspond to a sharp implementation of 
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the Lorentz gauge fixing, where one can simply set £[uv] = and suppress the corresponding rows 
and columns in the Hessian. 

Next we perform the TT decomposition on the symmetric part of the vielbein fluctuation 



d- 



(41) 



-V 2 



Ay) (~ v2 )- 



and the associated field redefinitions £ M — > y — V 2 — ^ £^ and cr — 

With these definitions, and dropping bars from the background quantities for notational sim- 
plicity, the quadratic part of the action (|36p is 



,(2) 
h T h T 

r (2) 
r (2) 



-v 2 + 



d(d-3) + 4 ^d-2 



-v 2 + 



d(d- 1) 
a(d- 2) - 1 



2d 



fl-(2-C)A 



d 



(d- 1) 2(d- 1) -a(d-2) 



2d 



ad 

v (d-2)(2- C )a-4 E 



ad (2 - C) 



,(2) 



,(2) 
ha 



2(d- 1) -a(d-2 
d - 2 2 (d 2 - 3d + 2) a - (d/3 - 2) 2 



-A 



a 



(42) 
(43) 
(44) 

(45) 
h 

w-iy (46) 

We notice that for ft = \ (a(d — 2) + 2) we can get rid of the mixed term. In the rest of the paper 
we will work in this "diagonal" gauge. In this case the trace part reduces to 



4d 



v 2 + 



da(d - 2) 

a(d-2) (d-4 + C) 



2(d- 1) - a(d- 2) 

da (d - 2 + C) 



2 (d 2 - 3d + 2) a - (d/3 - 2) 



;R-2 



-(d-2)a + d/3-2 (d-1) 



da 



d 



Vgh(-v 



R 



2 (d 2 - 3d + 2) a - (d/3 - 2) 



rA 



,(2) 



l d-2 
2 2d 



-V 2 + 



d-4 + C 



-R 



2d 



1 + 



c 



d-2 



A 



h 



2(d - 1) - a(d - 2) 2(d - 1) - a(d - 2) 
After decomposing the diffeomorphism ghost in its transverse and longitudinal parts, and ab- 
sorbing v 7 — V 2 in the latter, the ghost action is the sum of 



(47) 



The Lorentz ghosts do not propagate and following standard perturbative procedure one could 
neglect them entirely, but we will follow and introduce a cutoff for them too. The corresponding 
contribution to the FRGE is computed, together with other traces, in Appendix II. 



B. Beta functions 



The FRGE can now be calculated by introducing a cutoff separately in each spin sector. (This is 
known as a "type b" cutoff.) Using the same heat kernel methods that we have employed in section 
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IIII Al the expression for OfT^ can be expanded to linear order in R. Comparing the terms of order 
zero and one in R in the FRGE yields: 

/ 9 A \ k d 

- a '(iib)=i^ (B ' + " a) (49) 

where r\ = —dtG/G and Ai, B{ are in general polynomials in A = k~ 2 A. From here one can find 
the beta functions of the dimensionless parameters G = k d ~ 2 G and A: 

dt G=(d-2)G + -^- (50) 

8 t A = - 2A + gA 1 + 2B 1 A + G(A 1 B 2 -A 2 B 1 ) 

2(1 + B 2 G) 

In the following two sections we will give explicit results using specific cutoffs. 

For numerical results in d = 4 we will always use the optimized cutoff ([3]). For a discussion of 

n 

the dependence on the shape of the function Rk(z) we refer to [34j. We will instead concentrate on 
the differences between cutoffs of type I vs. II and type a vs b. For the type la case we refer again 
to the extensive discussion in [34], whose results we have checked independently. We will report in 
detail the results for the cases lb and lib, and highlight the differences with the case la. 



C. Type lb cutoff 



First we choose as reference operator, in each spin sector, the "bare" Laplacian —V 2 . The cutoff 
is a function Rk(— V 2 ). This is called a cutoff of type lb. The calculation of the coefficients A\, 
A 2 , B\ and B 2 for arbitrary dimension and cutoff shape is described in Appendix A. Here we just 
report the result in d = 4 and for the cutoff profile ([3]) : 
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The result is still quite general: it depends on the parameter £, which allows us to interpolate 
continuously between the purely metric formulation (£ = 0) and the purely tetrad formulation 
(C = 1)> on the arbitrary gauge parameter a, which allows us to test the gauge dependence of the 
results, and on the parameter /i that allows us to weigh differently the contribution of the Lorentz 
ghosts. 

Let us now describe the main features of these flows. Both in the metric and in the tetrad 
formulations, a UV-attractive fixed point is found for all values of \i and for a not too large. Its 
location and the corresponding critical exponents i? (which are defined as minus the eigenvalues of 
the stability matrix) are given in table I in the metric (( = 0) or tetrad (£ = 1) formalism, in the 
gauges a = 0, a = 1 and with two different values of the dimensionless Lorentz ghost parameter 
Ji = fj,/k: /I = oo, and \x = 1.2. The former corresponds to neglecting the Lorentz ghosts entirely 
and the latter is chosen to ease comparison with [34], who found that this value gives results that 
are closest to the metric formalism in the gauge and scheme they use. Note that the case £ = 0, 
a = 1 correspond to the purely metric flow with type lb cutoff, which had already been discussed 
previously in the literature. Indeed the second row in table I agrees with the third row in table II 
in fl. 

Whereas with a type la cutoff the fixed point becomes UV-repulsive, and a limit cycle develops, 
for fi sufficiently large, with the type lb cutoff it remains UV-attractive for arbitrarily large \x. 
This is a nice feature of this scheme, because it means that the fixed point can also be found if one 
adopts the perturbative prescription of neglecting the Lorentz ghosts entirely. However, the results 
in the tetrad formalism match most closely those of the metric formalism when fi is chosen to be a 
bit larger than one. As with the type la cutoff, for \x smaller than a critical value /2 C , the critical 
exponents become real. We find jl c ~ 0.705 for a = and jj, c ~ 0.715 for a = 1. 

The dependence of the universal quantities on the gauge parameters is illustrated in Fig. 1 . The 
slow decrease of the real part of the critical exponent for < a < 2 is in agreement with earlier 
calculations in the metric formalism (see e.g. fig. 9 in [7j|). The results of different schemes seem 
to converge for a — > which, we recall, is believed to give the physically most reliable picture. On 
the other hand when a is greater than some value of order 2 the fixed point becomes repulsive, 



reproducing the behavior that had been observed in [34] for large /L It is tempting to conjecture 
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Scheme 


A* 


G* 




<d 


lb, c = 


0, a = 







0.1569 


0.9028 


0.1416 


2.147± 2.620i 


lb, c = 


0, a ~ 


l 




0.1715 


0.7012 


0.1203 


1.689 ±2.486i 


lb, c = 


i, A = 


oo, a 


- 


0.2288 


1.363 


0.3119 


2.086 ± 2.042i 


lb, c = 


i, A = 


oo, a 


= 1 


0.2478 


0.9472 


0.2347 


0.595±3.753i 


lb, c = 


i, A = 


1.2, a 


= 


0.0691 


1.518 


0.1050 


2.237 ± 1.248i 


lb, c = 


l, A = 


1.2, a 


= 1 


0.0798 


1.3196 


0.1053 


1.892 ± 1.093* 



TABLE I: The non trivial fixed point in the type lb cutoff in metric (( = 0) and tetrad (£ = 1) formalism, 
in the gauges a = and a — 1 and with different weights of the Lorentz ghosts. Recall Re($)>0 implies 
that the fixed point is UV attractive. 
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FIG. 1: Plot of universal quantities as functions of the gauge parameter a for type lb cutoff. In the left 
panel the product A*G*, in the right panel the real part of the critical exponent. 

that also in the cutoff scheme la used in [34 1 the fixed point would have the usual properties, even 
for large fl, if one could choose a closer to zero. The strong /independence that had been observed 
there may be due to a particularly strong a-dependence. Altogether it appears that with a type 
lb cutoff, the tetrad formalism leads to results that are qualitatively similar to those of the metric 
formalism, and that the correspondence is best when < a < 1 and the Lorentz ghosts are turned 
on, with a parameter A that is a little larger than one. 



D. Type lib cutoff 



We call type lib a cutoff imposed separately on each spin-component of the graviton and taking 
as reference operator the Laplace-type operator that appears in the corresponding Hessian, includ- 
ing the curvature terms, but not the term proportional to the cosmological constant. The rationale 
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Scheme 


A* 


G* 




d 


lib, c = 


0, a = 









0.1052 


0.7216 


0.0759 


2.562 ± 1.566i 


lib, c = 


0,a = 


1 






0.0924 


0.5557 


0.0513 


2.424 ± 1.270i 


lib, c = 


l,ji = 


CO, 


a = 





0.1406 


1.0176 


0.1431 


2.595 ± 1.131* 


lib, c = 


1, fi = 


00, 


a 


1 


0.1369 


0.8427 


0.1154 


2.300 ±0.99H 


lib, c = 


1, p, = 


1.2, 


i \ = 





0.0394 


1.0008 


0.0398 


2.640 ±0.730i 


lib, c = 


l,fi = 


1.2, 


a 


1 


0.0361 


0.8299 


0.0300 


2.547±0.634i 



TABLE II: The non trivial fixed point in the type lib cutoff in metric (£ = 0) and tetrad (( = 1) formalism, 
in the gauges a = and a = 1 and with different weights of the Lorentz ghosts. 



for excluding the cosmological constant term is that the cosmological constant is a running coupling 
and if one included it in the reference operator, it would not remain fixed in the course of the flow. 
Here we choose a reference operator that remains fixed along the flow. 4 As we have already seen 
in the case of the fermions, the type II cutoffs tend to give somewhat simpler final formulae than 
the corresponding type I cutoffs, because to leading order one always finds traces of the function 
dtRk/Pk an d it is not necessary to expand the denominators. 



The coefficients A\, A2 are the same as with a type lb cutoff and are given in (j52|53j) . The 
coefficients B\ and B2 for arbitrary dimension and cutoff shape are given in Appendix III. In d = 4 
and for the cutoff profile Q, they become 
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2a|^A 
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(56) 
(57) 



! — C)A l-a(2-C)A 

Table II gives the UV-attractive fixed point and the corresponding critical exponents in the 
metric (£ = 0) or tetrad (C = 1) formalism, in the gauges a = 0, 1, and with two different values of 
the Lorentz ghost parameter, (l = 00, and (l = 1.2. 

The results are qualitatively similar to the ones obtained with the type lb cutoff. This is in 
line with all the results obtained previously in the Einstein-Hilbert truncation. The non trivial FP 
exists and has complex critical exponents for all values of (t greater than a critical value /2 C , which 
is approximately equal to 0.766 for a = and 0.748 for a = 1. For small (i the FP moves towards 
negative values of A. For large fx the fixed point remains UV attractive, in contrast to the result 



4 Cutoffs that depend on the full Hessian, including the terms proportional to the cosmological constant, were called 
of type III in [5|, where they have been applied to the Einstein-Hilbert truncation. 
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a 



FIG. 2: Plot of universal quantities as functions of the gauge parameter a for type lib cutoff. In the left 
panel the product A*G*, in the right panel the real part of the critical exponent. 

found in [34] with the type la cutoff scheme (which we have independently verified) . In particular 
we find that the FP has properties close to the standard ones of the metric formulation also when 
the Lorentz ghosts are neglected. Figure 2 gives the gauge-dependence of the universal quantities 
AG and Note that the real part of the scaling exponent $ is particularly stable in this scheme, 
for < a < 1. 



E. Type Ha cutoff 



For completeness we mention here also the results for the cutoff of type Ha, which had been 
discussed first in section IVC of [5]. In this scheme only the gauge a = 1 is easily computable. 
In this gauge it is enough to split the metric fluctuation into its trace and tracefree parts to write 
the Hessian of the Einstein-Hilbert action as two minimal Laplace-type operators. The cutoff 
is then defined as a function of these operators, including the curvature terms but excluding the 
cosmological constant term. This prescription leads to particularly simple expressions. 

The coefficients Ay , A<i are the same as with a other cutoff types considered here and are given 
in (I52|53p . The coefficients By and B<i in d = 4 and for the cutoff profile ([3]) are simply 

2 - 3C 27(2 - C) 



By 

B 2 



1 

127T 
1 

48^ 



l-(2 + C)A l-(2-C)A 
2 - 3C 27(2 - C) 



40 



+ B 1 



(58) 
(59) 



.l-(2 + C)A l-(2-C)A. 
These expressions coincide with (|56l57p when one puts a = 1 there. As a consequence, all properties 
of the flow are the same and we will not discuss this case further. 

It is nevertheless interesting to understand the reason for this coincidence, which is not restricted 
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to d = 4 and is also independent of the shape of the function R k . 5 For the sake of simplicity we 
shall discuss here only the case £ = 0, but the result is general. Since in all cases the trace field 
h is treated separately, and its contribution is the same for type a and b cutoffs, it is enough to 
consider the tracefree part of the graviton, namely the components hz£, £^ and a. In the gauge 
a = 1, the Hessian in the tracefree subsector is a minimal second order operator of the form 

- V 2 + C T R - 2A (60) 

where Ct = d ^(d-^ • When one uses a cutoff of type Ha (no further decomposition) the contribu- 
tion of this sector to the r.h.s. of the FRGE is 

E d t R(X n ) - r]R k {X n ) 
„ P fc (A n )-2A [ ' 

where A n are the eigenvalues of the operator O = — V 2 + CtR- One can divide these eigenvalues 
into three classes, depending on the spin of the corresponding eigenfunction. Upon using the TT 
decomposition (|4"T|) one finds that the eigenvalues of O on fields of type hj^ , V^£„ — and 
V^V^c — ^V 2 er are equal to the eigenvalues of the operators in square brackets in (I42I43I44|) . 
stripped of the A terms. We denote these operators O tt = —V 2 + CtR, = —V 2 + C^R and 
O u = —V 2 + C a R and the corresponding eigenvalues X^ T , A| and A^. So, the trace (|63l) is equal 
to 

d t R(Xl T )- V R k (Xl T ) ^ d t R(Xj)- V R k (Xj) dtRjX^-rjRkjK) 
2? P k (XD-2A + ^ P k (xi)-2A P k (K)~2A 1 ) 

Since for a = 1 the coefficients of A in (|42|43|44|) are all the same and equal to —2, this is recognized 
as the contribution of the fields hj^, £ M and a to the r.h.s. of the FRGE when one uses a cutoff of 
type lib. By a similar reasoning one also concludes that the ghost contribution is the same in the 
Ila and lib schemes. 

Things do not work in the same way for type I cutoffs, i.e. when the cutoff is a function of —V 2 . 
For a type la cutoff the contribution of the tracefree sector to the r.h.s. of the FRGE is 

E d t R{X n ) - r]Rk{X n ) . , 

t P k (X n ) + C T R-2K ' 1 ' 

where A n now denote the eigenvalues of —V 2 . This can be expanded as 

E d t R{X n ) - r]Rk(X n ) d t R(X n ) - rjR k (X n ) -A d t R(X n ) - i]R k (X n ) 

P k (X n ) + C T R-2K + P k (X n ) + C T R-2K + P k (X n ) + C T R-2K ' 1 ' 



5 The agreement between cutoffs Ila and lib for a = 1 had been noticed before in [38| . 
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where 22 > Z~2 anc ^ S denote the sum over eigenvalues of — V on hi„ , — and 

V^V^cr — ^V 2 cr respectively. On the other hand for a type lb cutoff the same contribution is 

TT £ 

E d t R(X n ) - r]R k (\ n ) d t R(\ n ) - r]R k {\ n ) d t R(X n ) - yRk(K) 

„ P fc (A„) + C T ii-2A ^P fe (A n ) + C^-2A ^ P *( A «) + ^ " 2A 

One clearly sees that the two traces are different. 



V. DISCUSSION 

The implementation of the FRGE in the presence of fermions and gravity presents some subtleties 
that had not been fully appreciated until recently. The sign ambiguity of the fermionic contribution 
to the running of Newton's constant had been known for a while, but it was regarded as just another 
aspect of the scheme dependence that is intrinsic to applications of the FRGE, albeit a particularly 
worrying one. Although a completely satisfactory understanding can only come from a treatment 
of physical observables, we have argued here that the correct treatment of fermion fields, when the 
Dirac operator is squared, is to use a cutoff that depends on — V 2 + 4r (type II cutoff). There 
also follows from our discussion that using a cutoff that depends on —V 2 (type I cutoff) may give 
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22j) have used 



physically incorrect results. Unfortunately, several earlier studies (in particular 
this scheme, so some of those results may have to be revised. We plan to return to this point in a 
future publication. 

Another issue is the use of tetrad vs. metric degrees of freedom. We have extended the analysis 



of tetrad gravity initiated in 34| by using a different cutoff (type lb and lib vs. type la) which 
allowed us to keep the diffeomorphism gauge parameter a arbitrary. We have found that the results 
for the running couplings using the tetrad formalism are qualitatively similar to those of the metric 
formalism, with some quirks. The following points should be noted. 

(i) For (" = one recovers the standard metric formalism. In the case of the type lb cutoff the 
results agree with the ones obtained earlier in the literature HQ, Q. The type lib cutoff with 
generic a had never been used before and the results obtained here are new. We have shown that 
for a = 1 they coincide with the ones found in [5| for the type Ha cutoff. For other values of a they 
differ only marginally from those obtained with other cutoff types and confirm the stability of the 
fixed point in the metric formalism. Type II cutoffs have the attractive feature that they lead to 
somewhat more compact expressions for the beta functions. 

ii) The case £ = 1 corresponds to the tetrad formalism. In this case a new ambiguity appears in 
the definition of the ghost sector: it can be parameterized by a mass \x that appears in the mixing 
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between diffeomorphism and Lorentz ghosts, or by the corresponding dimensionless parameter 
jl = n/k. This parameter is a priori arbitrary, but in order not to introduce additional mass scales 
into the problem it is natural to assume that it is of order one. On the other hand we recall that 
in perturbation theory and in the chosen gauge the Lorentz ghosts are neglected, since they do not 
propagate. This corresponds to taking jl = oo. 

iii) If one uses a type la cutoff and completely neglects the Lorentz ghosts, there is no attractive FP 
for positive G 34j. Instead, one finds a UV-repulsive fixed point surrounded by a UV-attractive 
limit cycle 6 . This is not the case when one imposes the cutoff separately on each spin component, 
as we did here. We find that both with type lb and lib cutoff's an attractive FP with complex 
critical exponents is present also when Lorentz ghosts are neglected, both for a = 1 and a = 0. 
This is reassuring because in the metric formalism the fixed point can be found even using the 
perturbative one loop beta functions. 

iv) If the contribution of Lorentz ghosts is added, as advocated in its effect is weighted by the 
parameter jl: it is strong for small jj, and weak for large jl. Since the ghosts are fermions, the fixed 
point is shifted towards negative A for decreasing jl. In addition, they have a systematic effect on 
the critical exponents: the modulus of the imaginary part decreases with decreasing jl and there is 
a critical value jl c under which the critical exponents become real. In the gauge a = 1, jl c = 0.715 
for a type lb cutoff, jl c = 0.748 for a type lib cutoff, and jl c ~ 0.8 for a type la cutoff. Similar 
behaviour is observed also for a = 0. On the other hand, the real part of the critical exponent 
decreases with increasing jl. With the type la cutoff this effect is most dramatic: the real part 
becomes negative for jl ~ 1.4, and this marks the appearance of the limit cycle. With the type 
II cutoffs discussed here the effect is much weaker and the fixed point becomes only slightly less 
attractive even in the limit jl — > oo, both for a = 1 and a = 0. 

v) The closest match between the tetrad and metric results is typically obtained if one chooses 
some value of jl that is not too far from one. This is always the case for the product AG, and in 
most cases also for the critical exponent. For the la cutoff this value was found to be approximately 
1.2. For the lb and lib cutoffs it is somewhat larger, depending on the quantities one is comparing. 
An exception occurs for the critical exponent in the case of a cutoff of type lib in the gauge a = 0, 
for which the best match occurs for jl — > oo. 

vi) Using type b cutoffs (i.e. decomposing the fields into irreducible components) has the advantage 
that one can keep the diffeomorphism gauge parameter a arbitrary. The gauge dependence of the 



This is different from the behaviour of the limit cycle discussed in [38 
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critical exponents is similar to what had been observed previously in the metric formalism, as long 
as a is not too much greater than one. In particular the real part of the critical exponent tends 
to decrease as a increases, for small values of a. In the limit a — )■ the /i-dependence becomes 
very weak and the critical exponents nicely converges towards a common value. On the other hand 
for a somewhat larger than one the fixed point either becomes complex or repulsive. This is the 
behaviour that had been observed in 



34| with the type la cutoff. This suggests that if we were able 
to compute the beta functions for this cutoff type with a/lwe would find that also with large \x 
the fixed point is present and has the standard properties for a sufficiently close to zero, 
vii) Altogether the results are very similar to those found in the metric formalism, except for the 
dependence on the new parameter /i, which is particularly strong for the type la cutoff. As argued 



m 



34j | . one can probably attribute the increased sensitivity of the results to the fact that in the 



tetrad formalism one has to deal with more unphysical degrees of freedom. The type la cutoff seems 
to be particularly sensitive to the off-shell, unphysical features of the flow. Type b cutoffs, where 
each spin component is treated separately, are less sensitive. 

viii) All of the preceding discussion is in the context of a "single metric truncation", i.e. one 



33|, the application of the 



assumes that the VEV of the fluctuation field is zero. As discussed in 
FRGE to gravity requires that the effective average action be considered in general a function of 
two variables. We plan to consider these more general truncations in a future paper. 

In conclusion, let us comment on the use of tetrad vs. metric variables. Since fermions exists 
in nature, it may seem in principle inevitable that gravity has to be described by tetrads. This 
would complicate the theory significantly. Every diffeomorphism-invariant functional of the metric 
can be viewed as a diffeomorphism and local Lorentz-invariant functional of the tetrad, but in the 
bimetric formalism there are many functionals of two tetrads that cannot be viewed as functionals 
of two metrics. Therefore, as already noted in [34j , the tetrad theory space is much bigger than the 
metric theory space. 

The necessity of using tetrads should, however, not be taken as a foregone conclusion. First of 
all, it is possible that the fermions occurring in nature are Kahler fermions. This would completely 
remove the argument for the use of tetrads, even in principle. Whether this is the case or not is a 
difficult issue that should be answered experimentally. For the time being one might just consider 
the use of Kahler fermions as a computational trick. Second, even if we stick to spinor matter, 
by squaring the Dirac operator and using a type II cutoff one can calculate many quantum effects 
due to fermions without ever having to use tetrad fields. The additional complications due to the 
Lorentz gauge fixing and the increased sensitivity to gauge and scheme choice advise against the 
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use of the tetrad formalism, as a matter of practical convenience. 
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VI. APPENDIX I: DIRAC SPECTRAL SUMS 



To compute the sum (J28J) we can use the Euler- Maclaurin formula 

p 



= I" F(x)dx-B r {F(n)+F(0))+ir (i^ 2 *- 1 ) („) - F^-^O) ^remainder (66) 

where Bi are the Bernoulli numbers. After collecting a volume contribution, the only terms we 
need to compute are the 0-th and 1-st power of R. Note that only the integral depends on R, and 
therefore, in dimensions d > 2 for the terms that we are interested in it is enough to compute the 
integral. 

Since the volume of the d— sphere is V(d) = ^jT + l) (47r)^ 2 ^ ^ H^ ) we on ^ have to 
isolate the terms in the integral proportional to R~ d l 2 and R l ~ d l 2 
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changing variables n — > n' — d/2 
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the terms we are interested in come from the integral of the two highest order power of n' 
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we can rewrite the sum Ylk=i \1 ~ ^) = li^i^ — 1) (^ — 2), and perform the integral 
dtRk 
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(68) 
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Collecting the volume of S d we obtain 
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d t Rk 
Pk 



- jl r d V(d) I k d - -k d ~ 2 R + O (R 2 ) 

r(| + i) (47r )f 1 ; V 24 ^ 



(71) 



23 



VII. APPENDIX II: TYPE IB CALCULATION 



We report here the detailed computation of the A and B coefficients of (I48I49P for a Type lb 
cutoff. The FRGE is the sum of traces over the irreducible components of the metric fluctuation 
defined in (l4~Tj) . They give: 
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Here a prime or a double prime indicate that the first or the first and the second eigenvalues 
have to be omitted from the^race (because ^ and a obey to some differential constraints, for 
more details see for example [5(). The contribution of the transverse and longitudinal parts of the 
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diffeomorphism ghosts are 
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Collecting the coefficients of J ^fg and — f y/gR we extract the A and B coefficients: 
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Here we have defined the dimensionless versions of the Q-functionals: Qd = k d Q± and Qd, = 

2 2 2 

2 

Finally let us consider the contribution of Lorentz ghosts. They do not propagate and therefore 
are usually neglected in the evaluation of the effective action in perturbation theory. Nevertheless 



if, following 
an amount 



34| . we impose a cutoff on their determinant they contribute to the r.h.s. of the FRGE 
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(82) 



Here we have introduced the arbitrary mass parameter [i (denoted in j^). In particular note 
that in the limit fj, — >■ oo the ghost contribution vanishes and one recovers the standard perturbative 
result where the Lorentz ghosts are neglected. Let A L and B L be the contribution of the Lorentz 
ghosts to the coefficients A\ and B±, defined in (|48l49p . From the above we read off 
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Note the appearance of Rk instead of Pk in the denominators. In general the Q functionals 
Q n ( fl^^/v^ ) Can ^ e com P u t e< ^ explicitly, with cutoff ([3]), in terms of hypergeometric functions. 
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For the calculations in four dimensions we only need the following 



VIII. APPENDIX III: TYPE IIB CALCULATION 



We report here the A and B coefficients of (|48l49p for a Type lib cutoff. The contributions of 
the irreducible components of the metric fluctuation to the r.h.s. of the FRGE are 
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The contribution of Lorentz ghosts is the same as in the Type lb case. The coefficients A± and A2 
are the same as in equations ([781 \7U\i , whereas 
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